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1.  INTRODUCTION 

The  funds  received  from  the  AFOSR  were  to  be  used  to  1)  purchase  materials  and  machine 
friction  test  specimens  and  2)  refurbish  a  Brown  and  Sharpe  Coordinate  Measurement  Machine 
(CMM).  The  specimens  were  used  in  an  initially  Carnegie  Mellon  University  (CMU)  supported 
and  then  US  Navy  supported  research  project  titled,  “Physical  Models  Of  Friction  Constraint 
And  The  Effect  Of  Friction  On  Bladed  Disk  Dynamics.”  For  the  first  year,  the  graduate  student 
working  on  the  project  was  supported  by  CMU  as  a  Teaching  Assistant.  The  AFOSR  funds 
allowed  us  to  make  specimens  and  begin  testing  before  we  started  getting  Navy  funding.  The 
CMM  machine  allowed  us  to  check  the  dimensions  of  the  specimens  that  we  had  manufactured 
and  characterize  how  wear  affected  the  contact  geometry.  Additional,  details  are  provided 
below. 


2.  TEST  SPECIMENS  AND  EQUIPMENT  PURCHASED 
2.1  Test  Specimens 

Approximately  $15,000  of  the  AFOSR  funds  were  used  to  manufacture  test  specimens  used  to 
measure  cyclic  response.  The  geometry  of  the  test  specimens  were  similar  to  that  of  under¬ 
platform  friction  blade  dampers  used  in  the  hot  section  of  gas  turbines.  We  made  specimens  with 
two  types  of  geometries.  One  type  had  spherical  heads  and  the  other  flat,  cylindrical  heads 
(Figure  1).  The  damper  specimens  fit  inside  “fork  specimens,”  an  example  of  which  is  shown  in 
Figure  2.  The  specimens  were  then  put  into  an  MTS  machine.  The  fork  was  attached  to  the  top 
actuator  in  the  MTS  machine  and  the  damper  specimen  to  the  bottom  actuator.  A  clamp  was 
then  applied  across  the  fork  that  produced  a  compressive  load  on  the  outside  of  the  fork  that  was 
transmitted  to  the  damper  specimen  (see  Figure  3).  Thus,  the  damper  and  fork  specimens  were 
connected  through  a  friction  interface. 


Damper  Specimen-Type  I  :  Cylindrical  Contact  Surface 


Damper  Specimen-Type  II :  Spherical  Contact  Surface 


Figure  1  Photographs  of  Internal  Damper  Specimens:  Top  photos  -  Flat  Cylinder,  Bottom 

photos  -  Spherical 
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Figure  2  Photograph  of  Outside,  Fork  Specimen 
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— Clamping  forks 
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|  Normal  Load 
Damper  head 


Actuator 

Figure  3  Left  Photo  Shows  Clamp  and  MTS  Machine.  Right  Schematic  Explains  the 

Arrangement  of  the  Test  Pieces 


The  normal  load  on  the  contact  surfaces  was  measured  using  a  strain  gage  mounted  at  the  center 
of  the  damper  specimen.  A  sinusoidally  varying  force  was  applied  in  the  vertical  direction  to  the 
frictionally  constrained  system  that  caused  the  friction  interface  between  the  fork  and  the  damper 
to  slip.  The  MTS  test  system  plotted  the  resulting  displacement  as  a  function  of  the  applied  load. 
Typical  sets  of  hysteresis  curves  are  shown  in  Figure  4  for  a  range  of  amplitudes  of  input  forces. 
The  figures  on  the  left  are  for  a  normal  load  (from  the  clamp)  of  2000  lbf  and  the  ones  on  the 
right  for  a  normal  load  of  4000  lbf. 
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Figure  4  Hysteresis  Curves:  (a)  and  (b)  are  for  the  flat,  cylindrical  specimen,  (c)  and  (d)  for  the 

spherical  specimen 

Tests  were  conducted  with  the  specimens  over  a  range  of  normal  loads  and  amplitudes.  The  data 
was  processed  and  compared  with  theoretically  predicted  results.  The  agreement  between  the 
theory  and  the  experimental  results  was  good.  A  detail  discussion  of  the  testing  and  analysis  was 
presented  as  a  paper  at  the  ASME  Turbo  Expo  2004  Power  for  Land,  Sea,  and  Air  in  Vienna, 
Austria  in  June  2004.  A  copy  of  the  paper  is  attached  as  an  Appendix.  The  paper  has  also  been 
accepted  for  publication  in  the  ASME  Journal  of  Engineering  for  Gas  Turbines  and  Power. 

2.2  CMM  Machine 

AFOSR  funds  were  also  used  to  refurbish  a  Brown  and  Sharpe  CMM  machine,  Figure  4.  The 
CMM  machine  was  used  throughout  the  project  to  check  the  dimensions  of  test  specimens  as 
initially  manufactured  and  also  to  characterize  geometry  changes  caused  by  wear.  We  received 
approval  from  the  AFOSR  project  engineer,  Dr.  D.  Mook,  to  use  the  funds  for  the  purpose. 


Figure  4  Refurbished  CMM  Machine 
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3.  FOLLOW-ON  RESEARCH 


The  quasi-static  testing  and  analysis  that  we  did  using  the  AFOSR  test  specimens  showed  that 
friction  contact  could  be  predicted  from  first  principles  using  nonlinear  contact  mechanics.  We 
have  subsequently  perfonncd  vibration  tests  on  turbine  blade-like  specimens  to  show  that  the 
quasi-static  results  apply  to  vibrating  under-platform  dampers.  Photographs  of  the  test  apparatus 
and  the  under-platform  damper  (with  spherical  heads)  are  shown  in  Figure  5. 


Figure  5  Vibration  Test  Apparatus  and  Under-Platform  Damper  Specimen 


The  performance  of  the  damper  is  depicted  in  Figure  6.  In  this  figure,  the  peak  vibratory 
amplitude  of  the  blade  tip  is  plotted  as  a  function  of  the  magnitude  of  the  excitation  force.  The 
curve  is  nonlinear  since  friction  effects  are  nonlinear.  The  experimental  data  agree  very  well 
with  analytical  predictions  indicated  by  the  continuous  curve.  The  analytical  prediction  were 
based  on  the  theory  developed  using  the  quasi-static  test  results  discussed  in  Section  2. 
Consequently,  this  follow-on  research  has  shown  that  damper  performance  can,  in  fact,  be 
predicted  from  first  principles  using  nonlinear  contact  mechanics. 

The  research  of  the  student  who  perfonned  these  studies,  Mr.  Kyun  Hee  Koh,  is  currently  being 
documented  in  his  PhD  theses.  The  more  recent  results  on  the  vibration  tests  were  just 
completed  this  summer  and  were  presented  by  Mr.  Koh  at  the  GUIde  Consortium  annual 
meeting,  July  2004,  in  Pittsburgh,  PA. 
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Figure  6  Damper  Performance  Curve  Compares  Prediction  with  Measured  Data 
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ABSTRACT 

This  paper  discusses  an  approach  for  characterizing  the 
dynamic  behavior  of  a  friction  damper.  To  accomplish  this, 
the  deflection  of  the  damper  is  measured  as  a  function  of  an 
applied  force  for  a  range  of  amplitudes,  normal  loads,  and 
excitation  frequencies.  The  resulting  hysteresis  curves  are  used 
to  generate  curves  of  nonlinear  stiffhess  and  damping  as  a 
function  of  the  amplitude  of  motion.  A  method  of  presenting 
this  information  in  a  dimensionless  format  is  demonstrated. 
This  format  allows  direct  comparisons  of  the  nonlinear  stiffness 
and  damping  of  actual  dampers  with  that  often  used  in 
analytical  models  to  compute  the  dynamic  response  of 
frictionally  damped  turbine  blades.  It  is  shown  that  for  the 
case  of  a  damper  with  a  spherical  head  significant  differences 
exist  between  the  actual  behavior  of  the  damper  and  that  often 
assumed  in  simple  analytical  models. 

In  addition,  Mindlin’s  analysis  of  a  sphere  on  a  half  space 
is  used  to  estimate  the  damper’s  stiffness  as  well  as  its 
theoretical  hysteresis  curves.  The  hysteresis  curves  are  then 
used  to  determine  dimensionless  stiffness  and  damping  curves. 
The  results  compare  favorably  with  those  found 
experimentally. 


INTRODUCTION 

Friction  dampers  are  used  in  gas  turbine  engines  to 
increase  energy  dissipation  and  reduce  vibratory  response. 
Dampers  are  pieces  of  metal  that  fit  between  the  blades  and  that 
are  held  in  place  by  centrifugal  force.  They  dissipate  energy 
when  sliding  occurs  across  the  friction  interface.  In  addition, 
the  damper  acts  like  a  weak  nonlinear  spring  that  affects  the 
frequency  of  peak,  vibratory  response.  Thus,  friction  dampers 
control  the  amplitude  and,  less  strongly,  the  frequency  of  the 
peak  response. 


^On  leave  from  Politecnico  di  Torino,  Turin,  Italy. 


Two  recent  review  papers  on  the  general  topic  of  friction 
damping  are  by  Feeny  et  al.  [1]  and  by  Gaul  and  Nitsche  [2], 
Predicting  and  optimizing  the  performance  of  turbine  blade 
friction  dampers  is  an  active  area  of  research  in  its  own  right. 
For  examples  consider  the  review  paper  by  Griffin  [3],  the 
papers  by  Yang  and  Menq  [4,  5,  6,  7],  Sanliturk  et  al.  [8], 
Petrov  and  Ewins  [9,  10],  Szwedowicz  et  al.  [11],  Cha  and 
Sinha  [12,  13],  Csaba  [14],  and  Sextro  et  al.  [15].  Computer 
codes  that  developed  from  these  research  efforts  are  widely 
used  by  the  gas  turbine  industry  to  predict  and  optimize  damper 
performance.  Nearly  all  of  the  codes  use  spring  elements  that 
can  slip  as  their  basic  friction  element  [16].  For  example  a 
friction  damper  that  connects  adjoining  blade  platforms  may  be 
modeled  by  a  spring  connection  (or  a  series  of  spring 
connections  in  parallel)  that  slips  if  the  magnitude  of  its  force 
equals  a  certain  limiting  value,  pN.  Analyzing  the  effect  of 
the  friction  force  becomes  more  complex  when  the  damper  and 
blade’s  contact  surfaces  are  inclined  so  that  the  normal  load 
varies  dynamically  during  a  cycle  [17,  18]  or  if  the  relative 
motion  across  the  contacting  surfaces  follows  a  two 
dimensional  path  [19,  20]. 

A  significant  effort  has  been  made  by  engine 
manufacturers  to  assess  how  well  the  resulting  computer  codes 
predict  damper  behavior.  Recently,  the  US  Air  Force  Office 
of  Scientific  Research,  the  US  Navy,  and  the  Sandia  National 
Laboratories  sponsored  two  workshops  on  contact  mechanics 
and  friction  damping  [21].  It  was  attended  by  representatives 
from  General  Electric,  Pratt  &  Whitney,  Rolls  Royce,  and 
Siemens  Westinghouse  as  well  as  representatives  from 
academia  and  US  government.  The  general  opinion  expressed 
by  the  industrial  representatives  at  the  workshop  was  that  the 
codes  used  by  industry  have  to  be  calibrated  with  test  data 
before  they  can  be  used  to  reliably  predict  the  performance  of  a 
particular  blade/damper  design  under  engine  conditions. 
Specifically,  the  consensus  view  was  that  the  damper  stiffness 
used  in  the  computer  code  that  correlates  with  the  test  data  is 
usually  significantly  less  than  that  calculated  from  a  finite 
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element  model  of  the  damper.  The  purpose  of  the  research 
reported  in  this  paper  is  to  measure  the  properties  of  a  real,  but 
simple,  damper  and  to  assess  why  it  might  be  difficult  to 
characterize  its  behavior  with  a  simple  spring  element. 

In  order  to  perform  this  assessment,  a  prototype  damper 
with  a  spherical  head  was  manufactured  and  tested  under 
various  conditions.  The  data  from  the  resulting  friction  force  - 
displacement  hysteresis  curves  were  processed  to  determine  the 
damper’s  effective  nonlinear  damping  and  stiffness  as  a 
function  of  the  amplitude  of  the  motion  across  the  friction 
interface.  The  results  were  then  compared  with  that  of  the 
simple  spring  damper  model  and  were  seen  to  exhibit 
significant  differences. 

The  contact  problem  of  a  sphere  on  a  half  space  subjected 
to  periodically  applied  shear  loads  was  analyzed  by  Mindlin 
[22,  23]  and  by  Deresiewicz  [24].  Experimental  studies  that 
support  the  theory  have  been  reported  by  Mindlin  [25],  Johnson 
[26],  and  Goodman  and  Brown  [27].  The  hysteresis  curves  of 
damper  with  a  spherical  head  in  contact  with  a  flat  plate  can  be 
predicted  analytically  using  Mindlin’s  analysis.  The  resulting 
curves  can  then  be  processed  in  the  same  manner  as  the 
experimental  data  to  predict  the  damper’s  effective  nonlinear 
damping  and  stiffness.  It  will  be  shown  that  Mindlin’s 
solution  accurately  captures  the  experimental  behavior  of  the 
damper  over  a  range  of  testing  conditions. 

The  organization  of  this  paper  is  as  follows.  First,  the 
analytical  properties  of  the  spring  damper  element  are 
reviewed.  They  will  be  used  as  the  benchmark  in  this  study. 
Of  specific  interest  is  the  manner  in  which  the  dimensionless, 
nonlinear  stiffness  and  damping  curves  are  used  to  characterize 
its  behavior.  Section  3,  will  describe  how  the  hysteresis  curves 
of  the  test  dampers  are  measured  and  then  processed  to  get 
dimensionless,  nonlinear  stiffness  and  damping  curves  that  can 
be  compared  with  those  of  the  spring  damper.  In  Section  4, 
Mindlin’s  solution  will  be  used  to  generate  theoretical  values  of 
nonlinear  stiffness  and  damping  that  are  then  compared  with 
the  experimental  results.  Implications  of  the  results  are 
discussed  in  Section  5.  Concluding  remarks  are  provided  in 
the  last  section. 


1.  SPRING  DAMPER  MODEL:  BENCHMARK 

To  understand  how  friction  dampers  work,  consider  the 
case  of  a  single  degree  of  freedom  oscillator  constrained  by  the 
spring  friction  damper  shown  in  Figure  1.  For  a  spring 
damper  the  energy  dissipated  during  one  cycle  of  motion  is  the 
area  enclosed  by  the  bilinear  hysteresis  curve  shown  in  Figure 

2.  It  is  considered  a  “macroslip”  element  since  it  does  not 
dissipate  any  energy  unless  there  is  gross  slip  across  the 
interface.  It  is  characterized  by  its  stiffness,  kj,  and  the  slip 
load,  N. 

If  the  nonlinear  friction  force  from  the  spring  damper  that 
resists  the  motion  is  fN ,  the  equation  of  motion  is 

m x  +  cx+kx  =  f0coscot- fN(t)  (1) 

The  nonlinear  friction  force  may  be  approximated  by  a 
spring  and  viscous  damper  acting  in  parallel,  i.e.. 


x 


Relative 


K 


Figure  1:  Spring  damper  element 


Figure  2:  Hysteresis  from  spring  damper  element 


fN(t)~Kx  +  cex 


(2) 


The  effective  stiffness  and  damping  are  found  using  the 
Harmonic  Balance  (HB)  method.  Some  examples  of  the 
application  of  HB  method  can  be  found  in  [16,  28,  29,  30].  In 
the  HB  method  it  is  assumed  that  the  displacement  is  sinusoidal 
and  that  the  nonlinear  force  has  the  same  period  as  the  motion. 
The  force  is  expanded  in  a  Fourier  series  and  truncated  after  its 
fundamental  terms  since  the  mass  spring  system  is  an  oscillator 
that  filters  out  the  higher  frequency  components.  The 
fundamental  terms  in  the  Fourier  series  can  be  expressed  so  that 
one  is  in-phase  with  the  motion  and  one  is  out-of-phase.  The 
in-phase  term  effectively  represents  a  spring  force  and  the  out- 
of-phase  term  an  equivalent  damping  force.  Since  the  Fourier 
coefficients  depend  on  the  amplitude  of  the  motion,  the 
effective  stiffness  and  damping  are  nonlinear  functions  of  the 
amplitude.  The  amplitude  of  the  motion  can  be  normalized  by 
djviding  by  the  distance  required  to  cause  slip,  i.e., 
A  =  A!{iiN / kd)  .  Then  the  non-dimensionalized  effective 
stiffness  and  damping  of  the  spring  friction  element  are  given 
by  the  equations [31]: 


m  i 


n 


cos 


n 


0.5sin<2cos  1  1- 


cM)  = 


ce(A) 

kdlnct) 


(3) 

(4) 


where  the  arc  cosine  is  defined  to  lie  between  0  and  n .  The 
dimensionless  stiffness  and  damping  of  a  spring  friction 
element  are  plotted  as  a  function  of  amplitude  in  Figures  3  and 
4.  They  completely  characterize  the  dynamic  behavior  of  the 
spring  damper  and  determine  the  response  of  the  single  degree 
of  freedom  system.  Note,  that  for  A  <  1 ,  the  damper  does  not 
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2.  EXPERIMENT 


slip.  Consequently,  the  damper’s  effective  stiffness  remains 
equal  to  its  elastic  value  and  its  effective  damping  is  zero,  i.e., 
the  spring  element  exhibits  only  macroslip  behavior. 


Figure  4:  Nondimensionalized  effective  damping 
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Figure  5:  Experimental  setup  for  quasi-static  test 


2.1  Experimental  setup 

Two  different  experiments  are  used  to  measure  the 
behavior  of  the  dampers.  One  is  a  quasi-static  test  conducted  at 
0.1  Hertz.  The  second  is  a  dynamic  test  conducted  at  50 
Hertz.  They  are  described  in  the  following  sub-sections. 

Quasi-static  test.  The  quasi-static  tests  were  done  using 
the  MTS  testing  machine  shown  in  Figure  5.  The  specimens 
have  two  parts.  A  fork  type  specimen  with  two  prongs  is 
mounted  in  the  upper  part  of  the  MTS  and  damper  type 
specimen.  Figure  6,  in  the  lower  part.  The  friction  damper 
specimen  is  inserted  between  the  prongs  which  are  then 
compressed  with  a  clamp.  The  prongs  of  fork  are  thin  plates 
that  are  meant  to  represent  the  blade’s  platform.  Strain  gages 
are  attached  to  the  flat  faces  of  the  damper  specimen.  They 
are  used  to  monitor  the  normal  load  on  the  friction  interface 
during  the  tests.  An  extensometer  is  used  to  measure  the 
relative  displacement.  A  special  fixture  was  designed  so  that 
the  extensometer  measures  the  relative  displacement  directly 
across  the  friction  interface.  This  eliminates  errors  caused  by 
the  rotational  deflection  of  the  fork  and  improves  the 
repeatability  of  the  measurement.  The  load  cell  in  the  MTS 
machine  is  used  to  measure  the  friction  force. 

In  order  to  facilitate  testing,  the  specimens  were  made 
larger  than  those  commonly  used  in  engines.  If  the 
dimensions  of  a  damper  are  scaled  by  a  factor  S  then  the 
displacements  are  also  scaled  by  S  and  are  easier  to  measure. 
To  compensate,  the  normal  loads  are  scaled  in  proportion  to  the 
square  of  S  in  order  to  maintain  realistic  contact  stresses.  The 
specimens  are  made  of  Maraging  steel.  The  advantages  of 
Maraging  steel  are  its  fair  machinability  before  heat  treating 
and  its  high  yield  stress  after.  Once  it  is  heat-treated,  it  has  a 
yield  stress  of  2,068-2,758  MPa,  comparable  to  that  of  the 
super  Nickel  alloys  used  in  turbines.  An  advantage  of  using 
this  high  strength  alloy  is  that  the  deformation  in  the  contact 
region  tends  to  stay  elastic.  The  faces  of  the  damper  that  come 
in  contact  with  the  prongs  of  the  fork  are  spherical  and  have  a 
radius  of  curvature  of  127  mm. 

Dynamic  test.  The  apparatus  used  in  the  dynamic  test  is 
shown  in  Figure  7.  Its  complete  description  is  given  in  [32]. 
In  this  experiment,  the  horizontal  beam  element  is  excited  by  a 
shaker  at  50  Hertz.  A  friction  damper  is  attached  to  the  beam 
at  its  center.  The  damper  consists  of  a  hemispherical  damper 
specimen  in  contact  with  a  flat  plate.  Its  radius  of  curvature  is 
12.7  mm,  an  order  of  magnitude  smaller  that  those  tested  in  the 
MTS  machine.  The  normal  load  is  applied  using  pulleys  and  a 
dead  weight.  The  normal  load  was  maintained  at  69  N  during 
the  tests  reported  here.  The  damper  shear  force  is  measured 
by  a  load  cell  in  series  with  the  damper  specimen.  The 
relative  motion  across  the  friction  interface  is  measured  with  a 
laser  vibrometer. 

Hysteresis  curves  were  generated  under  dynamic 
conditions  using  this  apparatus.  The  tests  provide  an 
interesting  contrast  to  the  quasi-static  experiment  since  the 
frequency  is  two  orders  of  magnitude  higher,  the  radius  of 
curvature  one  order  of  magnitude  lower,  and  the  normal  load 
two  orders  of  magnitude  lower  than  in  the  quasi-static  tests. 
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Figure  6:  Friction  damper  specimen 
for  quasi-static  test 


Laser  vibrometer 


N 
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Figure  7:  Experimental  setup  for  dynamic  test 
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Figure  8:  Hysteresis  curves  from  quasi-static  tests 


2.2  Measurement 

The  tests  were  conducted  for  a  range  of  amplitudes  and 
at  different  normal  loads.  Representative  results  for  quasi¬ 
static  test  are  shown  in  Figure  8.  The  quasi-static  tests  were 
done  for  two  different  normal  loads,  8.9  kN  and  17.8  kN. 
Similar  curves  were  generated  using  the  dynamic  test 
apparatus,  Figure  9. 

2.3  Data  processing 

The  hysteresis  curves  from  the  tests  were  processed  in  a 
manner  consistent  with  that  used  to  determine  the  equivalent 
stiffness  and  damping  of  the  spring  damper.  First,  it  was 
necessary  to  non-dimensionalize  the  amplitude  of  motion. 
This  requires  knowing  the  damper’s  physical  stiffness  kd  and 
the  force  required  for  gross  slip,  /rN  .  The  stiffness  was 
inferred  from  the  hysteresis  curves  by  averaging  the  slopes 
during  the  initial  unloading  and  loading  portions  of  the  cycle 
kd  ~  (A'j  +  k2)fl ,  Figure  10.  The  gross  slip  load  /^N  was 
also  determined  from  the  hysteresis  curves  by  applying 
sufficiently  large  amplitude  that  the  curves  flattened  out.  Thus, 
the  characteristic  length  A0  =  //TV  /  kd  was  determined  from  the 
experimental  data.  The  dimensionless  amplitude  A  could 
then  be  calculated  by  dividing  the  physical  amplitude  by  A0 . 

The  dimensionless  effective  stiffness  and  damping  are 
given  by: 

=  f  f„<x*Od0  (5) 

nAfiN 

^A)  =  -A—  [’fNsmOd0  (6) 

AfiN  * 

These  quantities  were  calculated  using  the  force  and 
displacement  information  from  the  hysteresis  curves  and 
integrating  the  data  numerically  in  MATLAB.  The  integration 
was  fairly  accurate  since  a  large  number  of  data  points,  about 
200,  were  measured  each  cycle. 

2.4  Key  results 

Figure  11  shows  the  nonlinear  effective  stiffness  and 
damping  plots  produced  by  the  quasi-static  and  dynamic  tests. 
Observe  that  the  quasi-static  experiments  with  the  two  different 
normal  loads  8.9  kN  and  17.8  kN  give  almost  identical  curves. 
The  normalized  data  points  from  the  dynamic  test  also  look 
quite  similar.  This  is  quite  an  interesting  result  since  it  tends 
to  imply  that  a  unique  set  of  dimensionless  stiffness  and 
damping  curves  exist  for  a  particular  damper  geometry  even 
though  the  test  conditions  varied  significantly  from  one  test  to 
the  next. 

When  our  experimentally  determined  curves  of  effective 
stiffness  are  compared  with  those  of  the  spring  damper  they 
agree  reasonably  well  for  large  displacements,  A  >  3  . 
However,  because  of  microslip  the  dimensionless  effective 
stiffness  of  the  experimental  damper  is  smaller  than  that  of  the 
spring  damper  when  the  amplitude  is  small.  More 
importantly,  the  experimental  data  has  much  more  damping  at 
low  amplitudes  than  the  spring  model  and  somewhat  less  at 
high  amplitudes. 
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Figure  9:  Hysteresis  curves  from  dynamic  tests 


Figure  10:  Calculating  damper  stiffness 
from  the  hysteresis 


Table  1:  Stiffness  values  from  experiments 


Dynamic  condition 
(Hz) 

Radius 

(m) 

Normal  load 
(N) 

kE.p(N/m) 

Dynamic  (50) 

0.0127 

69 

31.6x10s  :»** 

Quasi-static  (0.1) 

0.127 

8900 

327x10s 

Quasi-static  (0. 1 ) 

0.127 

17800 

401x10s 

In  the  next  section  a  classical  elasticity  solution  will  be 
investigated  in  order  to  assess  the  extent  to  which  it  captures 
the  behavior  seen  in  the  tests. 


3.  AN  ANALYTICAL  MODEL: 

MINDLIN’S  APPROXIMATE  SOLUTION 

In  this  section  an  approximate  solution  by  Mindlin  [25]  is 
summarized  for  the  elasticity  problem  of  a  tangential  load 
applied  to  a  sphere  in  contact  with  a  half  space.  It  will  be  used 
to  calculate  the  damper’s  stiffness,  kd ,  the  hysteresis  curves, 
and  the  dimensionless  effective  stiffness  and  damping  curves 
for  a  sphere  on  a  half  space. 

An  interesting  aspect  of  three-dimensional-contact 
problems  in  elasticity  is  that  the  relative  displacement  between 
the  contacting  bodies  approaches  a  finite  limit  in  the  far  field. 
For  example,  when  a  sphere  contacts  the  flat  surface  shown  in 
Figure  12,  this  behavior  limits  the  relative  displacement 
between  the  bodies.  Thus,  when  shear  loads  are  applied  the 
displacement  of  the  center  line  has  the  behavior  indicated  by  x 
in  Figure  13  provided  there  is  no  gross  slip. 

Consequently,  there  is  a  unique  stiffness  kx  that 
characterizes  the  far  field  behavior.  From  Mindlin ’s  theory  the 
far  field  stiffness  is: 


where  G  is  the  shear  modulus,  v  is  Poisson’s  ratio,  R  is  the 
radius  of  curvature,  a  is  the  contact  radius,  and  N  is  the  normal 
load. 

Mindlin  also  gives  analytical  expressions  that  define  the 
various  parts  of  the  hysteresis  curves,  Figure  14(a). 
Specifically,  the  relative  motion  is  given  by: 


A 


Figure  11:  Experimentally  generated  curves  compared 
with  benchmark 


The  values  of  damper  stiffness  determined  from  the 
experiments  are  summarized  in  Table  1.  Note  that  the 
stiffness  of  the  quasi-static  specimen  increased  with  normal 
load. 
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These  formulae  were  used  to  generate  the  hysteresis  curves 
shown  in  Figure  14(b).  The  curves  were  processed  in  the  same 
manner  as  the  experimental  data  to  determine  the  dimensionless 
effective  stiffness  and  damping  curves  for  the  Mindlin  model. 

N 


i 


Figure  12:  Friction  contact  of  sphere  on  a  flat  surface 


Figure  13:  Relative  displacement  between  two 
bodies  in  contact  ( F  <  pN  ) 


Figure  14:  Hysteresis  from  Mindlin  analysis 


Figure  15:  Comparison  between  experiments 
and  Mindlin  analysis 


Comparisons  between  the  experimental  results  and  the  result 
from  Mindlin’s  analysis  are  shown  in  Figure  15.  In  general, 
the  agreement  is  quite  good. 

Mindlin’s  estimate  of  the  damper’s  physical  stiffness  (7)  is 
a  function  of  material  properties,  normal  load,  and  radius.  A 
comparison  of  the  stiffnesses  measured  in  the  experiment  and 
those  calculated  from  (7)  is  shown  in  Table  2.  It  is  clear  that 
Mindlin’s  formula  certainly  captures  the  trend  seen  in  the 
experiments.  This  is  quite  remarkable  considering  the  wide 
difference  in  the  test  conditions  in  the  quasi-static  and  dynamic 
tests. 


Table  2:  Stiffness  comparison  between  experiments 
_ _ and  Mindlin  theory _ _ _ 


Dynamic 

Condition 

(Hz) 

Radius 

(m) 

Normal 
load  (N) 

kElp(N/m) 

k  Mindlin 

(N/m) 

Dynamic  (50) 

0.0127 

69 

31.6x10*  :"i* 

31.7x10* 

Quasi-static 

(0.1) 

0.127 

8900 

327x10* 

345x10* 

Quasi-static 

(01) 

0.127 

17800 

401x10* 

435x10* 

4.  MODEL  IMPLICATIONS 

Since  the  Mindlin  analysis  agreed  well  with  the 
experimental  data,  it  was  used  to  examine  other  implications  of 
the  theory.  First,  rewrite  (9)  -  (11)  in  dimensionless  form 
where  x  =  x/A0  and  fN  s  fN/(iN .  Then 


3  r  /  -  \?/3' 

—  J  for  OPQ 


Equation  (12)  implies 


fN’=/;{A)=iM-^A 

\ 


for  OPQ 

(12) 

x  =  A  and 

$•* 

ii 

2  -V/3 

1A 

(13) 
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for  STP  (15) 


determine  how  the  stiffness  and  damping  change  for  small 
motions.  To  first  order,  the  dimensionless,  nonlinear  stiffness 
decreases  by  A! 6  and  the  damping  equals  2Af9  .  In 
dimensional  terms,  for  example,  this  means  that  for  physical 
amplitudes  A  less  than  jiNlkdX\\c  effective  damping  is 
approximately  equal  to 


c.(A)  = 


2  Ah] 
9a)7t/uN 


(24) 


Inverting  (14)  and  (15)  to  find  explicit  expressions  for  the 
cyclic  force  fN  results  in 

fN{X  *)=/»' -2  [|^cos^+(i-/;f3+i 31  J(16) 

for  PRS 

M*  ^)=-/;+2{i-(0  -^cos^+(i-/;f+i *  j  (17) 

for  STP 

where  fN*  =  fN*\A]  is  given  by  (13).  Now  substitute  (16) 
and  (17)  into  (5)  and  (6).  It  is  clear  that 


~  f2  fN(A,  6\  cos  OdO 

kM)=^ — - 

A 

(18) 

% 

c 

‘w 

cL 

* 

ii 

/—N 

> - " 

(19) 

Since  the  0  dependence  in  (18)  and  (19)  is  integrated  out, 
it  is  found  that  the  dimensionless  effective  stiffness  and 
damping  are  only  functions  of  dimensionless  amplitude  A  . 

-  -  fc(a) 

A 

(20) 

caa: >= — 

A 

(21) 

where 

0)  COS  Odd 

(22) 

F'W=“f  A (4  0)sin0d0 

(23) 

Consequently,  it  is  clear  that  the  dimensionless  effective 
stiffness  and  damping  predicted  by  Mindlin’s  theory  are 
independent  of  size,  material,  normal  load,  and  friction 
coefficient.  Thus,  according  to  Mindlin’s  theory,  there  are 
unique  nonlinear,  dimensionless  stiffness  and  damping  curves 
for  a  sphere  on  a  half  space. 

Another  result  from  Mindlin’s  theory  is  that  you  can 


5.  CONCLUDING  REMARKS 

While  in  general  it  is  not  possible  to  use 
nondimensionalization  techniques  for  highly  non-linear 
systems,  in  the  case  of  the  spherical  damper  geometry 
considered  in  this  paper  it  was  found  that  the  test  data  and  the 
analytical  solution  indicate  there  is  a  unique  set  of 
dimensionless  curves  that  characterize  its  performance.  In 
addition,  two  additional  parameters,  kd  and  p,  need  to  be 
known  in  order  to  complete  the  description  of  a  specific 
damper.  It  was  found  that  the  stiffness  of  the  spherical 
damper  kd  could  be  reasonably  well  approximated  by 
Mindlin’s  theory,  equation  (7).  Since  Mindlin’s  theory  is 
based  on  an  elasticity  analysis  that  takes  into  account  nonlinear 
contact  this  results  suggests  that  the  stiffness  of  dampers  used 
by  industry  might  be  calculated  from  elasticity  theory  provided 
sufficient  care  is  taken  to  model  the  effect  of  nonlinear  contact 
and  the  resulting,  very  local  stress  concentrations. 

The  harmonic  balance  method  was  used  to  calculate  the 
vibratory  response  of  a  frictionally  damped,  single  degree  of 
freedom  oscillator  in  which  the  friction  element  was 
represented  by  a  spring  damper  [16]  and  then  by  the  Mindlin 
model.  The  results  are  shown  in  the  damper  optimization 
curves  of  Figure  17.  Damper  optimization  curves  plot  the 
maximum  amplitude  of  response,  Apeak,  as  a  function  of  the 
normal  load,  N,  applied  to  the  friction  joint  [16].  In  most 
applications,  turbine  blade  dampers  are  designed  to  work  in  the 
white-colored  range.  In  fact,  this  range  corresponds  to  the 
white  colored  range  in  Figure  18  where  the  dimensionless 
amplitude  A  varies  between  0  and  1 .094.  In  this  operating 
range,  it  is  clear  that  the  real  damping  was  very  different  from 
that  predicted  by  the  spring  model.  The  results  imply  that 
microslip  effects  may  be  very  important  over  the  part  of  the 
operating  range  typically  used  by  friction  damper  designers. 
Based  on  these  observations,  it  is  believed  that  the  problem  that 
the  friction  damper  designers  have  in  calculating  the  correct 
damper  stiffness  is  due  in  part  to  the  fact  that  their  models  may 
not  properly  account  for  microslip.  As  a  result,  they  are 
changing  a  parameter  in  their  macroslip  model,  1^,  in  order  to 
try  to  match  behavior  that  is  dominated  by  microslip  effects. 
This  hypothesis  that  industrial  dampers  primarily  work  in  the 
microslip  regime  is  also  supported  by  the  fact  that  industrial 
dampers  rarely  show  gross  wear  after  use  in  the  engine. 
However,  whenever  the  test  data  were  taken  for  A  >  2, 
significant  levels  of  wear  were  observed  after  only  a  few 
hundred  cycles. 

Lastly,  an  additional  challenge  remains  in  predicting  the 
coefficient  of  friction,  p.  During  the  experiments  it  was 
observed  that  p  changed  with  time  as  the  surface  finish  of  the 
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Figure  16:  Dimensionless  hysteresis  loop 
from  Mindlin  theory 


Figure  17:  Design  range  in  damper  optimization  curve 


Figure  18:  Design  range  in  effective  damping  curve 


contact  surfaces  degraded.2  Consequently,  it  is  an  open  issue 
what  its  value  would  be  after  the  damper  has  operated  for  some 
time  in  the  engine.  In  our  tests  \i  increased  by  about  a  factor 


Note  that  we  had  to  periodically  measure  p  and  take  its  change  into 
account  in  developing  the  dimensionless  damping  and  stiffness  curves  quoted 
in  this  paper. 


of  two  when  the  motion  across  the  friction  interface  became 
large  after  enough  cycles.  If  the  damper  does,  in  fact,  work  in 
the  microslip  regime,  then  from  (24)  it  is  clear  that  an  increase 
in  p  would  result  in  less  damping  and  higher  vibratory 
response.  Consequently,  it  is  important  to  insure  that  p 
achieves  a  stable  value  representative  of  engine  operating 
conditions  in  evaluating  the  damper’s  effectiveness  during  early 
testing. 
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